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^ : INTRODUCTION 

p ' The overwhelming majority of theoretical papers on phase transition of hadronic 

^ ■ matter H to quark gluon plasma ( H ^ QGP) ignores constituent quarks Q which 

we shall call below briefly valons (the term proposed by R.Hwa). The present paper is 
fundamentally based on the conception of those valons as real entities having the same 
^ \ quantum numbers as current quarks and the mass which was more than once calculated 

d \ theoretically as some 300 MeV and used for explanation of experiments (e.g., see^). The 

possibility of existence of a special state of strong interacting matter, that of deconfined 
valon gas with still broken chiral symmetry, is discussed here. 

Such a possibility, i.e. appearance of the third, Q phase intermediate between H 
and QGP had been discussed already in papers^"^ (we apologize for possibly missing 
some other papers). The idea may be traced back to E.Shuryak who pointed out that 
deconfinement and restoration of chiral symmetry might not coincide^. However in the 
papers mentioned above the result turned out pessimistic. E.g., on the /x — T diagram 
(/i - chemical potential, T - temperature) possible Q phase (allowing also for admixture 
of pions necessary first of all as Goldstone particles) occupies but a very small area at 
large T and rather small // (see e.g. Fig.l taken from^).The "constituent quark phase 
is strongly suppressed"^. 

Meanwhile the intuitive physical pattern^'^ based on the bag model ideology, al- 
though being extremely rough, points to the quite a different possibility schematically 
shown in Fig.2. In fact, a nucleus consists of nucleons which in turn are bags contain- 
ing valons, which are nothing else than a massive clouds of strong interacting current 
quarks g, antiquarks q and gluons g with quantum numbers of a single quark. If we 
compress the nucleus (Fig. 2a) until nucleons become close packed and fill up the entire 
space (Fig. 26) they merge into a single bag with valons deconfined (still with chiral 



symmetry broken) and freely moving within the common bag (Fig. 2c). Further com- 
pression leads to a new situation, that of close packed and merging valons (Fig.2(i). 
Now q, q and g which earlier were confined within valons can freely propagate within 
a common bag (Fig.2e). Here they are deconfined and massless, chiral symmetry is 
authomatically restored and this very state is QGP. 

Thus we expect two phase transitions to take place: the one, H Q, and (at 
higher density and/or temperature) another, Q ^ QGP. 

In this paper we are going to show, in distinction to^~^, that with reasonable 
values of bag parameters Bq for the valon phase and Bq for QGP there actually appears 
possibility for a pattern shown in Fig.2, i.e. the possibihty of the overall H <-> QGP 
transition to proceed only via the intermediate Q phase. 

The first (H ^ Q) transition establishes limitations on possibility of existence 
of normal hadrons. This is the reason for interpreting it as the Hagedorn limit and to 
consider the temperature of transition, i.e. of valon deconfinement, T^, as the Hagedorn 
temperature Th- The second transition restores chiral symmetry and is characterized 
by the temperature Tch- Thus the Hagedorn temperature appears to be not merely a 
very rough approximation to chiral symmetry restoration/ breaking temperature T^h as 
is frequently believed. It has its own physical meaning of temperature of constituent 
quark (=valon) deconfinement: Th = Tj^. 



THERMODYNAMICAL FORMALISM 

We are treating phase transitions following traditional manner fundamentals 
of which were given by Hagedorn^. As the first step we consider the partition func- 
tion Z'^{T, ii,V) depending on T, and system volume V for each phase in the first 
approximation (which means: for pointlike particles): 
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where gff '^and mj are degeneracy factors and masses for Bose and Fermi i-th type par- 
ticles respectively for each j-th phase (j designates hadronic i^, valonic Q, and current 
{QGP), q phases respectively). lnZ°(T, //,y) is a result of summation over particles 
which play considerable role (under T and /j in question) in the given phase. In partic- 
ular, for hadronic phase H we have taken into account nucleons (protons and neutrons, 
M = 940 MeV), pions (m^ = 140 MeV) and A hyperons (tba ^ 1150 MeV). The 
constituent phase Q contains light valons = constituent quarks (m„~m(i~320 MeV), 
some admixture of strange valon quarks (ms~512 MeV), and again pions. The term 
InZ^^g stands here for Q and QGP phases and reflects effective interaction of quarks 
and gluons with the QCD-vacuum. This interaction in Q and QGP phases is assumed 
as 
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where BQ,Bg are some constant bag parameters differing from each other (to be dis- 
cussed below). 

Then by usual differentiations of Zj{T, fi,V) we obtain for each phase energy 

density ef\ pressure pf\ and particle density nf^as functions of internal parameters 
T,fi,V. 

Having this as a basis and following the procedure proposed^ ^ ( more exactly 
in^ ) we introduce mutual interaction of particles by ascribing them hard core radii in 
the Van der Waals manner. Namely, we substitute with real particle density n'^^^ 
and then obtain "real" energy density e^^^ and pressure p*^''^ for each of three phases 
{j = H,Q,QGPy. 

where each of the parameters vj is hard core volume characterizing j-th phase particles: 



p-(r,/x)= — , (6) 



= ^vrr^ , (7) 

'"Q = l^^Q > (8) 

current quarks and gluons arc of course pointlikc (vg = 0). 

At given T and V of the system the most preferable (i.e. stable) phase is the one 
having the largest pressure pj and smallest /ij. At equilibrium of each two phases theire 
pressures and effective chemical potentials are to be equal. Three types of phase equi- 
librium curves are to be analyzed: the one corresponding to the valon deconfinement 
transition {H Q) 

I^H = 3//q; PniTd, i^h) = PqiT^, i^h/^) , (9) 
direct transition (H ^ QGP) 

Hh = 3//^; Ph{Tc, I^h) = Pq{Tc, I^h/^) , (10) 
and chiral transition (Q <-> QGP) 

I^Q = l^q; PqiTch, I^q) = PqiT^h, I^h/^) ■ (11) 

The point of coexistence of all three phases (i.e. the triple point) is defined by the 
following conditions: 

ph(t*, /.#) = pq(t*, ^J*) = p,(r#, /.#) (12) 

Before turning to study of equilibrium conditions and phase transitions let us discuss 
the chosen set of parameters. 



BASIC PARAMETERS 



There are two pairs of decisively important parameters: firstly, particle sizes, e.g. 
average root square radii of nucleons N,rN, understood as (r^Ar)^/^, and valons, rg, 
which designates (r^^)^/^. Secondly bag pressures, Bq for the Q phase and Bg for 
QGP. 

Let us discuss sizes. They are used^~^ to describe interaction of particles in a 
simplified form by considering them as black balls forming a Van der Waals gas. This 
means that we approximate their interaction potential f/(r) as f/ = 00 for r < and 
U = for r > Tc, Tc being the core radius (in other calculations we have used also some 
continuous potential U{r)). In^ the nucleon interaction was instead described by two 
versions of the continuous U{r) satisfactorily describing nuclear properties. It seems 
sufficient in this paper to consider and Q simply as having cores^'^. For nucleon 
radius and volume rjv, vn we usually assume, as most resonable, the values: 

vn ^ 0.8fm; vn ^ 2.13fm=^ (13) 

(however sometimes we shall consider also — 0.7 fm). Since the nucleus has the 
radius Ra = * A^^^ {A - atomic weight) with tq — 1.10 — 1.15fm, the fraction of the 
nucleus volume occupied by nucleons is 

(TN/rAf ~ 0.26 ^ 0.38 . (14) 

Thus it seems sufficient to compress a nucleus only 3-4 times to come to close packing 
of nucleons. 

The size of valons has been estimated^° from analysis of properties of various 
hadrons within the Additive Quark Model (AQM); when expressed by the square root 
radius rQ(for — 0.8 fm) this estimate sounds as: 

0.20 < tq < 0.36fm. (15) 

We shall concentrate on some middle resonable value, tq = O.Sfm, although other values 
within limitation (15) will be also used (let us remark that estimating tq within black 
ball valon model and AQM from QQ collision cross section ctqq = 7r(2rQ)^ ~ {l/9)aNN 
where is the A^'A'" collision cross section taken, at not too high energies, as 40 mb, 
we get rg = 0.2 fm. For tevatron energies it is much larger). 

Now, about Bq and Bg. The choice of their values is much more disputable, 
being at the same time very essential. Previous authors^"^ were (seemimgly) guided by 
the results of lattice calculations in which confinement/deconfinement (c/d) and chiral 
symmetry breaking/restoration (b/r) transitions had been found to coincide at = 0. 
Moreover "deconfinement" was (and still is) usually understood as the straightforward 
H QGP transition. Accordingly, e.g. a condition was superimposed^ on pressures 
Ph and Pg of H and QGP phases at the transition point (for = 0) : ph{Tc,0) = 
Pq{Tc, 0). Herefrom the ratio (3 — Bg/ Bq was obtained thus leaving only one of these 
two parameters undefined. In other works similar coincidence of {c/d) and (b/r) at /x = 
was also taken as a basis. Herefrom, for a freely chosen Bq, followed the value of Bg. 
All authors accordingly used mainly the values (3 = Bg/ Bq ~ 3-j-4 with Bq more or less 
close to the MIT bag value, Bmit ^ 56 MeV/fm^ ( Eg = 56 ^ 150 MeV/fm^ Bg = 
200 ^ 500 MeV/fia, keeping in all cases^-^ ^ ~ 3 ^ 4 ). 

Assumption concerning Bq ~ Bmit seems to be reasonable at least in the vicinity 
of the H ^ Q transition. In fact, here valon density in Q phase, at least in the beginning 
of the valon deconfinement process, is the same as in a nucleon (see Fig. 26). Accordingly 
we shall use the assumption: 

50 < Bg ~ Bmit < 100 MeV/iia. (16) 



On the contrary, argumentation which have lead to the above mentioned esti- 
mate of f3 = Bq/Bq ~ 3 is not admissible for us since difference between confine- 
ment /deconfinement (c/d) of valons (which takes place with broken chiral symmetry) 
and chiral symmetry breaking/restoration {b/r) is exactly the phenomenon we are look- 
ing for. Accordingly this difference should not be excluded at the start of the analysis 
as has been, as a matter of fact, done e.g. in^. Thus Bg is a parameter which needs its 
own physical foundation. Since it refers to QGP, it should be tightly connected with 
the vacuum pressure Byac — 500 -j- 1000 MeV/fm^ and be close to it. Accordingly we 
are going to use larger values of P — Bq/Bq putting it equal to /? ~ 10. Below we shall 
use the following set of reasonable parameters as the standard one: 

Bg = 500 MeV/fm^; Bq = 50 MeV/fm=^; r^v = 0.8fm; rg = 0.3fm. (17) 

The influence of Bq and Bg as well as rj variations are partially investigated. 

MAIN RESULTS 

Let us start with extreme cases: /i = 0, T ^ 0(Figs.3,4) and T = 0, fi ^ 0(Fig.5). 

Fig. 3 clearly demonstrates characteristic pattern of double phase transition at 
the p — T plane for chosen standard parameter set (17). It is evident that two phase 
transitions are well separated, Q is the most stable one (i.e. its pressure pg is higher 
than that of other phases) within the temperature interval of AT = T^h — ~ 50 
MeV, with absolute values of transition temperatures being of the order of the direct 
phase transition {H ^ QGP) temperature T^: ~ ISOMeV, Td ~ U5Mev, T^h ^ 
195Afet'. The value of this "temperature corridor" of Q phase depends on Bg (see 
Fig.4(6)) so that it becomes larger (up to 90 MeV) for larger values olBg (when Bg = 
1 GeV/fm^,/3 ~ 20, Tcu ~ 230 MeV) and shrinks to a triple point (which are marked 
by the symbol within all the figures) T* ~ 140 MeV for most frequently choosen 
Bg ~ 200 MeV/im /9 ~ 4, close to that in^ The same effect stays for Bg dependence: 
the closer it is to Bg, the narrower is T-region allowing for Q phase existence, with 
the biggest value (corresponding to the triple point T* ~ 180 MeV) Bg ~ 160 MeV 
(i.e. 13 ~ 4.5). For all combinations olB's the chiral symmetry restores at T^h which is 
always larger than the critical temperature of direct transition = Th^qgp obtained in 
calculations when ignoring existence of valons. Note that we do not show the influence 
of r parameters becouse it is negligeable in this T — /i region. 

The similar pattern is presented in Fig. 5 for T = 0, 7^ 0. for the standard 
parameter set(17). Here again there is a possibility for appearance of Q phase region 
but now it is greately suppressed: if /9 = 10 while Bg has its minimal reasonable value 
50 MeV/fm^ the covered by Q phase on {p, fi} diagramm is rather small and it 
shrinks into a triple point already at Bg ~ 80 MeV/fm^; this region becomes of course 
broader for larger (3 values. However it is essential that valons tell even at T = 0. For 
more details concerning T = see 

Now we go over to the general case of T, /i ^ 0. The {/i, T} diagramm for the 
standard parameter set (17) is shown in Fig. 6. It is clear that there exists a wide 
corridor around the curve for direct H ^ QGP transition which remains as well for 
various combinations of parameters {Bg, Bg,rN, and rg) (see Figs. 7). It shrinks and 
even disappears with // growth only at rather definite conditions, namely for small /3(the 
shape of Q phase region for /3 < 3 shown in Fig. 7(a), as well as in Fig.l, reminds a thin 
banana) and for extremely large rg (Fig. 7(6)). Otherwise transition from H phase to 
QGP proceeds only via valon phase (for detailes see ^^). 



Note that under not too large fj, < 1 GeV the master parameter is /3 = B^/Bq 
while influence of particle interaction (i.e. choice of values of r^- and tq) is quite 
negligeable. Symbol line in Fig. 7(6) corresponds to the case of all particles taken point 
like. It differs considerably from the real particle diagram only at ~ 1 GeV, and 
covers the area having the shape of some exotic fruit (" pineapple" ) . s 

Critical values e*^*" for energy densities ( that is the ones on phase equilibrium 
curves ) are presented in Figs. 8. Let us point out (see Fig.8(a)) that the latent heat, 
i.e. the jump Ae*^ = (e^ — ejj) at the deconfinement transition is much smaller than 
that at the chiral one, Aech = {<^f' — ), and much smaller than that for direct tran- 
sition ignoring valons (Fig.8(6)). The decomfinement transition requires rather "soft" 
experimental conditions. 

DISCUSSION AND SUMMARY 

Of course, obtained numerical values, in particular those of and Tch, should not 
be taken literary. E.g., they essentially depend on Bq and Bg values chosen in a rather 
rough manner. According to e.g. Shuryak for hadronic substance at T = 0, -Be// 
(efficient Bgff describing the summarized influence of interaction of particles between 
themselves and the vacuum) is not reducible to two constants but changes continuously 
increasing with baryon number density ub from some lower value of the order of Bmit 
up to -Be// ~ B^ac for QGP. However, our calculations show that baryon number 
density in Q phase, ng, is much closer to that in the H phase than to the one in the 
QGP. Thus the use of only two values - Bq and Bg - seems to be not too bad. 

Further, our calculations show that particle density of H phase in the vicinity of 
transition curve to Q phase at relatively large fi > 1 GeV is very close to its limiting 
value ~ (close packing) and use of the Van der Waals gas approximation for it is 
very suspicious. But the same weak point is met in the papers which lead to pessimistic 
conclusions concerning possibility of the Q-phase^"^. Repulsive potential for nucleons 
considered in ^ seems of course to be more realistic; similar potential will be discussed 
in our paper This does not change the final result qualitatively. 

Nevertheless qualitative and even semiquantitative results seem to deserve atten- 
tion. Comparing Figs. 6, 7 with Fig.l we sec that for sufficiently large values of Bq/ Bq 
the pattern declared in the beginning of this paper (Fig.2) is almost fully supported. 
We may conclude that there actually exists possibility of hadronic phase to QGP trans- 
fer (at seemingly reasonable choice of parameters) to proceed only via an intermediate 
valon, i.e. constituent quark plus pions, phase. Therefore, we expect existence of two 
phase transitions. For = transition temperatures differ by some 50 MeV. One of 
them, deconfinement of valons, proceeds at the Hagedorn temperature = Th, above 
which, in the deconfined valon phase, hadrons exist but as a small thermal admixture. 
The second phase transition at Tch ~ 200 MeV corresponds to deconfiment/confinement 
of current quarks and gluons and, simultaneously, restoration/breaking of chiral sym- 
metry. Thus difference between and Tch is physically meaningful. 

The Q phase with deconfined constituent quarks (valons) having mass ~ 320 
MeV, due to mg/T > 1, is nonrelativistic, its state equation is that of nonrelativistic 
Van der Waals gas (since tq is small). It is only in the vicinity of the phase transition 
to QGP that valon size begins to tell. However, this last transition at not too large 
f3 = Bq/ Bq takes place long before close packing of valons is attained^^. This result is 
rather interesting and has more physical sense than had been put into the model. In 
fact the model of hard black balls, being quite sufficient for nucleons due to their high 



stability, is hardly valid for constituent quarks which may become partialy transparent 
and larger. Overlapping of valon tails securing exchange of current quarks and gluons 
between valons may happen at distances exceeding the assumed black ball radius. Thus 
the transition to QGP phase does not necessarily require for close packing of valons. 
This is not so for nucleons at H ■(-^ Q transition where the transition happens long 
after the considerable close packing is attained; this may be ascribed to necessity of 
overcoming the surface tension of nucleons. However for more detailed analysis it is 
desirable to consider more realistic repulsive potential for interaction of valons than 
simple hard core model (e.g. at the same manner as it has been done in^^). 

The obvious first apphcation of these results of course is to be to the evolution 
of matter formed at collision of highly relativistic heavy nuclei. Let us stress two most 
important items. 

First, as it have been shown above (Figs. 8) the latent heat value for H ^ Q 
transition is much smaller than that for the "upper" transition Q QGP (and much 
smaller than people had obtained for direct, i.e. H <-> QGP transition). The same 
is true^^ for transition pressure and ns values. This means that first deconfinment 
transition may proceed under much more "soft" conditions than it was expected for 
direct transition H ^ QGP: it is sufficient to compress the nucleus only 3-4 times 
applying not very high pressure. Thus this transition may happen at not too high 
energies of colliding nuclei, i.e. even at Bevalac or in Dubna^'^ where appearance of 
pure QGP is hardly possible. Accordingly, e.g. in the Dubna experiments one may 
expect to observe some effects impossible for simple pp collisions. In particular, our 
estimates show that admixture of A particles may attain in central rapidity region 
some 10 percents of nucleons and this reminds us Okonov group Dubna experiment 
in which A's were especially studied. 

Secondly, transition pattern described above should influence essentially the char- 
acter of hydrodynamical expansion and cooling of initially hot matter both if it takes 
place either in Q or in QGP phase. For QGP initial state the following picture seems to 
be natural. The system experiences the phase transition to QGP phase. After coohng 
down to chiral transition temperature Tch ~ 200MeV, it would further expand and cool 
down within a large temperature interval according to nonrclativistic state equation. 
Since at present ( and later when RHIC will be operating ) the attainable initial tem- 
perature may exceed Tch but slightly, the time spent for getting Tch, as well as duration 
of cooling down in H phase from ~ 150MeV to freezing temperature, Tf ~ 140 
MeV, may be smaller than the time of existence of Q phase with its nonrelativistic 
state equation. Therefore the existence of Q phase should tell itself essentially. Here 
detailed calculations for resulting rapidity distribution, etc. are necessary to compare 
with experiment (at present they are in progress). 

The results obtained seem to contradict to the statement based at lattice calcula- 
tions that at II — two phase transitions are to coincide. The situation here is however 
disputable. Although existence of Q phase was seemmgly never especialy looked for 
in lattice calculations there are some works in which existence of "heavy modes" was 
assumed in order to get rid of some results looking physically unsatisfactory^® ( the 
QGP phase does not have properties of an ideal relativistic gas ). Lattice specialists 
should answer the question whether the double phase transition with Q mass at 320 
MeV could be noticed in lattice calculations performed until now and how it should 
show itself. 
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FIGURE CAPTION. 

Figure 1. Phase transition curves in the T—n plane(from^): Q-phase region is between 
dashed and dotted lines. 

Figure 2. Schematic picture of nuclear matter transformation under compression. 

Figure 3. Phase diagram for /i = case at the P — T plane for standard parameter 
set (17): the solid, dashed and dotted lines correspond to the the pressures of hadronic 
(H), valonic (Q) and current q{QGP) phases respectively. 

Figure 4. The same plot as in Fig.3 for standard parameter set (solid lines) and 
varying B's: (a) - Bq ~ 100 (dashed hue), 160 MeV/fm^ (dotted line); (b) -Eg ~ 1000 
(dashed line), 200 MeV/fm^ (dotted line);"#" marks the triple points. 

Figure 5. Phase diagram in the p — /jl plane for T=0 and standard parameterset; 
designations are the same as in Fig.3. 

Figure 6. Phase transition curves in the T — n plane: valon deconfinement (dashed 
line), direct transition ignoring valons (solid line), chiral transition (short dashed). 

Figure 7. Transition curves bounding the Q-phase region: (a) - for standard 
and Bq, Bq = 50 MeV/fm^ (solid hues), 100 MeV/fm^ (long-dashed hues) and 150 
MeV/fm^ (dotted lines); short-dashed line corresponds to i?^ = 1 GeV/fm^; (b) - for 
standard S's and varying rg — 0.2 fm (sohd hues), 0.4 fm (dashed hues) and Vj = (f 
lines) . 

Figure 8. Energy densities ej curves for: (a) H and Q matter at deconfinement 
transition (solid lines), Q and QGP matter at chiral transition (dashed lines); (b) H 
and QGP matter at direct transition ignoring valons (dotted lines). 
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